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Abstract 


A  simple  ray  tracing  approximation  for  the  high 
frequency  scattering  from  a  two-dimensional  cavity  is 
developed.  Whereas  many  other  cavity  scattering  algorithms 
are  very  time  consuming,  this  method  is  very  swift. 

The  analytical  development  of  the  ray  tracing  approach 
is  performed  in  great  detail,  and  it  is  shown  how  the  radar 
cross  section  (RCS)  depends  on  the  cavity's  length  and  width, 
and  the  radar  wave's  angle  of  incidence.  This  explains  why 
the  cavity's  RCS  oscillates  as  a  function  of  incident  angle. 
The  RCS  of  a  two-dimensional  cavity  was  measured 
experimentally,  and  these  results  were  compared  to  computer 
calculations  based  on  the  high  frequency  ray  tracing  theory. 
The  comparison  was  favorable  in  the  sense  that  angular  RCS 
minimums  and  maximums  were  predicted  exactly  even  though 
absolute  accuracy  of  the  RCS  magnitude  decreased  for  incident 
angles  far  off-axis.  Overall,  once  the  method  is  extended  to 
three  dimensions,  this  technique  shows  promise  as  a  fast 
first  approximation  of  high  frequency  cavity  scattering. 
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ftgh  frequency  estimation  of 

TWO-DIMENSIONAL  CAVITY  SCATTERING 

I .  Introduction 

In  recent  years,  the  United  States  Air  Force  has  devoted 
much  effort  to  low  observables  technology.  Popularized  in 
the  media  by  '-.he  term  "stealth,"  this  area  of  research  has 
created  a  revolution  in  electronic  warfare  strategy.  Whereas 
the  survivability  of  combat  aircraft  was  previously  thought 
to  depend  on  sophisticated  jamming  equipment,  the  emphasis 
has  now  shifted  to  the  design  of  aircraft  capable  of  flying 
undetected  through  enemy  radar  coverage.  However,  the  task 
of  reducing  an  aircraft's  radar  "echo"  is  far  from  trivial. 
Until  recently,  this  echo,  commonly  referred  to  as  radar 
cross  section  (RCS),  could  not  be  accurately  predicted;  it 
required  actual  measurement.  This  approach  is  woefully 
inadequate  for  the  design  of  modern  aircraft  because  the 
iterative  construction,  measurement,  and  modification  of  an 
airframe  wastes  both  time  and  money.  An  ideal  solution  would 
be  a  computer  program  capable  of  estimating  the  RCS  of  any 
given  shape.  This  would  provide  the  engineer  with  the 
additional  advantage  of  being  able  to  try  many  more 
approaches  to  a  particular  design  problem  before  settling  on 


a  final  solution. 


Unfortunately,  no  single  comprehensive  airframe  RCS 
calculation  algorithm  has  been  developed  because  the  task  of 
predicting  the  RCS  of  even  simple  structures  is  staggeringly 
complex.  For  now,  the  building  block  approach  is  being  used 
to  work  toward  the  eventual  completion  of  a  comprehensive  RCS 
calculation  project.  One  of  the  basic  airframe  structures 
which  is  critical  to  this  process  is  the  engine  inlet  cavity. 
It  can  dominate  the  aircraft's  RCS  from  several  different 
aspects,  and  the  calculation  of  its  RCS  is  not  trivial  for 
even  the  simplest  cavity  designs.  Electromagnetic  scattering 
from  cavities  has  been  extensively  studied.  Moll  and  Seecamp 
analyzed  an  inlet  modelled  by  a  open-ended  cylindrical 
waveguide.  The  field  entering  the  cavity's  opening  is 
represented  by  a  series  of  waves  traveling  inward  while  the 
backscattered  field  is  a  series  of  waves  moving  from  the 
cavity's  closed  end  toward  the  opening.  The  sum  of  these  two 
fields  is  then  equated  to  the  total  field  at  the  cavity's 
opening,  and  the  Stratton-Chu  integral  is  used  to  solve  for 
the  scattered  fiela  (4). 

Johnson  and  Moffatt  used  the  Wiener-Hopf  technique  to 
find  the  scattering  matrices  for  the  cavity's  open  end,  and 
the  Green's  function  and  Moment  Method  were  employed  to  find 
che  matrix  for  the  termination  of  the  cavity.  The  RCS  is 
found  after  calculating  a  backscatter  coefficient  for  the 
entire  structure  (2). 
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And  finally,  Cooper  used  the  method  of  coupled  azimuthal 
potentials  and  appropriate  Dirichlet  boundary  conditions  to 
solve  for  the  fields  across  the  cavity's  opening  (2). 

Although  these  are  only  several  of  the  studies  which 
analyze  the  cavity  scattering  problem,  they  encompass  most  of 
the  basic  techniques.  However,  these  methods  typically 
suffer  from  one  serious  drawback.  They  all  tend  to  be 
computationally  intensive.  A  simpler  and  faster  technique 
could  have  some  applications  even  if  it  was  not  as  accurate 
as  any  of  the  methods  described  above.  A  simple  approach 
could  even  be  used  as  a  fast  first  approximation  for  other 
more  time  consuming  techniques  when  great  accuracy  is 
required.  Furthermore,  the  majority  of  the  existing  cavity 
scattering  methods  analyze  cy 1 indrical ly  shaped  cavities, 
whereas  many  inlets  of  critical  concern  to  RCS  engineers  have 
rectangular  openings. 

The  ray  tracing  method  for  the  analysis  of  cavity 
scattering  described  in  this  thesis  is  extremely  fast,  and  it 
applies  primarily  to  cavities  with  rectangular  openings.  It 
is  limited  to  the  high  frequency  case  where  the  cavity 
opening  is  much  larger  than  the  wavelength  of  the  incident 
electromagnetic  energy.  Also,  in  an  effort  to  simplify  this 
preliminary  attempt  at  developing  ray  tracing  theory,  the 
analysis  is  limited  to  two-dimensional  cavities  (see  Chapter 
III).  Once  the  success  of  this  technique  has  been 
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II.  Analytical  Development  of  the 
Ray  Tracing  Approach 

Introduction 

Before  the  complex  scattering  mechanisms  present  in 
three-dimensional  cavities  can  be  understood,  it  is  both 
instructive  and  informative  to  develop  a  model  to  analyze  the 
scattering  from  a  two-dimensional  cavity.  Although 
understanding  the  cross  section  of  a  two-dimensional  cavity 
provides  no  direct  benefit  since  these  structures  are  rarely 
encountered,  the  analysis  techniques  are  simple  a  r.  d 
straightforward.  Once  an  accurate  and  easily  implemented 
method  for  predicting  the  RCS  of  the  simple  two-dimensional 
cavity  can  be  developed,  these  same  techniques  can  eventually 
be  expanded  to  treat  a  three-dimensional  case  as  well. 
Although  a  three-dimensional  theoretical  model  is  the 
ultimate  goal,  an  accurate  tv/o-dimensional  model  is  a 
necessary  and  important  step  toward  achieving  that  end. 

First,  a  definition  of  a  two-dimensional  cavity  is  in 
order.  As  shown  in  Fig.  2-1,  the  cavity  has  two  finite 
dimensions  (L  and  w)  but  it  is  infinite  in  the  third 
dimension  (perpendicular  to  the  page).  The  direction  of 
m.cident  electromagnetic  radiation  is  shown  as  9.  This  angle 
can  vary  from  -90  to  90  degrees  but  only  within  the  plane  of 
the  page.  This  restriction  will  limit  the  necessary 
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When  the  effective  length  exactly  equals  three-quarters 
of  a  cavity  wavelength,  all  of  the  incident  radar  energy  is 
backscattered.  This  phenomenon  is  shown  in  Fig.  2-12. 

K 

T 

W 

1 

t L q .  2-12.  Total  Backscatter 

All  the  incident  rays  enter  the  cavity  above  this 
boundary  ray,  so  all  the  energy  is  backscattered. 

Case  £ 

This  case  occurs  when  the  cavity's  effective  length  is 
between  three-quarters  and  one  full  cavity  wavelength  (3  * /4 
<  Le  <  x  ).  Fig.  2-14  shows  a  case  4  cavity  with  its 
boundary  ray.  All  incident  rays  entering  below  this  boundary 
ray  will  be  backscattered. 


backscattered 

bundle 


Fig.  2-13.  Case  4  Cavity 
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Case  3 


This  case  occurs  when  the  effective  length  is  between 
one-half  and  three-quarters  of  a  cavity  wavelength  (  \/2  <  Le 
<  3  »/4). 

Fig.  2-11  shews  a  Case  3  cavity  vith  its  boundary  ray. 
The  backscatter  bundle  enters  the  cavity  above  this  ray,  and 
once  again  the  percentage  of  the  total  power  that  is 
backscattered  is 


where 


%  backscatter  =  d 
w 


tan  ©  =  d/ 2 

Le  -  W2 

or 


d  =  2(Le  -  *  /  2 )  tan© 


so 


d  2  (Le  -  * / 2 ) tan© 

%  backscatter  =  _  =  _ 

w  w 
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is  reflected  in  the  -9  direction,  and  none  is  backscattered. 
This  is  illustrated  in  Fig.  2-10.  Any  ray  entering  above 
this  boundary  ray  will  be  reflected,  so  every  ray  entering 
the  cavity  will  exit  at  an  angle  of  -9.  As  previously  noted, 
this  instance  of  total  backscatter  occurs  when 

Le  =  */2 


and 


?r 


tan  9  =  w  =  w__ 
*/2  Le 


9 


arctan 

w 

=  arctan 

w 

*  /  2 

Le 

- 

Fig.  2-10.  Total  Reflection 


backscattered 

bundle 


Fig.  2-9.  Case  2  Cavity 


percentage  of  the  total  power  that  is  backscattered  is  given 


%  backscatter  =  d 


where 


tan  9  =  w  -  d/ 2 
L0 


therefore 


d  =  2 (w  -  L0tan  9) 


%  backscatter  =  21- 


L>etan  9 


Once  again  there  is  a  limiting  case  for  a  cavity  of  this 
length.  When  the  effective  cavity  length  exactly  equals  one- 
half  of  a  cavity  wavelength,  all  of  the  incident  radar  energy 
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and 


Fig.  2-8.  Total  Backscatter 


tan  0  =  w/ 2  =  w/ 2 
A/ 4  Le 


or 


9  =  arctan 

2w 

=  arctan 

w 

\ 

2  Le 

1 

It  should  be  noted  that  this  total  backscatter  occurs  only 
when  Le  exactly  equals  */4.  However,  as  Le  aproaches  */4 
the  percentage  of  backscatter  approaches  100%. 


Case  2 

This  case  occurs  when  the  cavity's  effective  length  is 
between  one-quarter  and  one-half  cavity  wavelength.  (  */4  < 
Lp  <  a. / 2 ) .  Fig.  2-9  shows  a  case  2  cavity  with  its  boundary 
ray. 


In  this  case,  any  incident  ray  entering  below  this 
boundary  ray  will  be  backscattered ,  and  once  again  the 


Fig.  2-8.  Total  Backscatter 


and 


tan  0  =  w/ 2  =  w/  2 
x/ 4  Le 


or 


0 


arctan 


* 

2w 

=  arctan 

w 

X 

2  Le 

It  should  be  noted  that  this  total  backscatter  occurs  only 
when  L0  exactly  equals  */4.  However,  as  Le  aproaches  x/4 
the  percentage  of  backscatter  approaches  100%. 


Case  2 

This  case  occurs  when  the  cavity's  effective  length  is 
between  one-quarter  and  one-half  cavity  wavelength.  (  x/4  < 
<  x / 2 ) .  Fig.  2-9  shows  a  case  2  cavity  with  its  boundary 
ray . 

In  this  case,  any  incident  ray  entering  below  this 
boundary  ray  will  be  backscattered,  and  once  again  the 
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the  entire  power  which  is  backscattered  is 

%  backscatter  =  d 
w 

It  can  also  be  seen  that 

tan  0  =  d/ 2 

Le 


therefore 


and 


d  =  2  Le  tan  0 


2  Le  tan  0 

%  backscatter  =  _ 

w 

So,  for  the  Case  1  cavity,  the  percentage  of  backscatter 
increases  as  0  increases. 

When  the  effective  cavity  length  exactly  equals  one- 
quarter  of  a  wavelength,  an  interesting  phenomenon  occurs: 
all  the  radar  energy  is  backscattered.  This  can  be  seen  in 
Fig.  2-8  which  shows  the  boundary  ray  for  this  case.  The 
boundary  ray  lies  exactly  at  the  lower  edge  of  the  cavity 
opening,  so  all  rays  entering  above  it  (all  rays  entering  the 
cavity)  are  backscattered.  This  instance  of  total 
backscatter  occurs  when 

Le  =  a/4 
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Fic.  2-6.  Reflection  and  Backscatter 


Fig.  2-7.  Boundary  Ray 

Whereas  the  entire  cavity  has  an  opening  of  width  w, 
that  portion  of  the  opening  which  admits  backscattering  rays 
has  width  d  as  shown  in  Fig.  2-7.  Therefore,  that  portion  of 
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The  maximum  radar  cross  section  of  the  cavity  occurs 
when  9=0  degrees.  The  entire  back  end  of  the  cavity  is  in 
direct  view  at  this  aspect.  It  is  well  known  that  the  radar 
cross  section  of  a  flat  plate  for  high  frequencies  at  normal 
incidence  is: 

o(9  =  0)  =  4  t  A2 


where 

A  =  area  of  the  flat  plate 

-  wavelength  of  the  radar  energy 

This  RCS  is  the  maximum  return  from  the  cavity,  and  it  occurs 
exactly  at  0  =  0.  As  the  incident  angle  is  varied,  this  peak 
RCS  is  modulated  by  two  separate  phenomena.  The  first  is 
caused  by  the  changing  effective  aperture  area,  and  the 
second  is  the  varying  tradeoff  between  the  two  possible  exit 
directions  (0  and  -9). 

As  the  incident  angle  9  increases,  less  radar  energy  is 
allowed  to  enter  the  cavity  causing  the  RCS  to  decrease. 
This  is  shown  in  Fig.  2-4.  The  effective  aperture  area  is 
proportional  to  the  cosine  of  the  incident  angle. 

As  already  noted,  the  radar  cross  section  of  a  flat 
plate  is  proportional  to  its  area  squared,  so: 

o  =  4TTAoff  =  4"TT(Acos9)^  =  4  TTA^cos^B 


therefore 


\  =  2w 
tan  0 


It  is  important  to  note  that  this  cavity  wavelength  is  a 
function  of  incident  angle.  Therefore,  a  new  cavity 
wavelength  must  be  calculated  for  each  incident  angle  of 
interest.  Then  the  maximum  integer  number  of  cavity 
wavelengths  must  be  subtracted  from  the  total  length  of  the 
cavity  to  yield  the  effective  cavity  length.  A  further 
example  of  this  process  is  displayed  in  Fig.  2-3.  The  cavity 
is  over  two  wavelengths  long  as  shown  in  Fig.  2-3 (a),  and 
once  these  two  wavelengths  are  subtracted  from  the  overall 
length,  the  effective  cavity  is  shown  in  Fig.  2-3 (b). 


Fig.  2-3. 


Effective  Cavity  Length 


Fig.  2-2.  Wavelength  of  a  Two-Dimensional  Cavity 

scattering  using  this  ray  tracing  analysis  technique. 
Following  one  ray  of  the  incident  energy,  it  can  be  seen  in 
both  diagrams  that  entry  and  exit  locations  and  directions 
are  identical. 

The  distance  of  a  complete  cycle  of  bounces  as  shown  in 
Fig.  2-2  (a)  is  labeled  as  a  "wavelength."  This  is  not  to  be 
confused  with  the  actual  wavelength  of  the  electromagnetic 
wave.  The  cavity  wavelength  is  very  important  in  the 
simplification  of  this  analysis.  It  can  be  seen  in  Fig.  2  - 
2(a)  that 

tan  0  =  W  =  2w 
a/z  a 

where 

a  =  cavity  wavelength 
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Fig.  2-1.  General  Two-Dimensional  Cavity 


algebraic  equations  to  a  two-dimensional  plane.  It  is  this 
limitation  which  make  the  analysis  of  the  two-dimensional 
case  considerably  less  complicated  than  the  three-dimensional 
cavity. 

It  can  be  shown  that  as  a  first  order  approximation  all 
energy  entering  the  cavity  from  angle  ©  will  exit  at  either 
that  same  angle  9  or  its  mirror  image  (-©).  The  ratio  of  the 
powers  exiting  in  the  two  directions  depends  on  the  cavity's 
width  and  length  as  well  as  the  incident  angle  0. 

But  first,  the  structure  of  a  two-dimensional  cavity  can 
be  simplified  somewhat.  As  shown  in  Fig.  2-2  (a)  for  certain 
incident  angles,  the  entering  radar  energy  may  "bounce" 
several  times  within  the  cavity  before  exiting.  If  this 
energy  bounces  off  each  side  wall  before  it  hits  the  cavity 
end,  then  the  cavity  can  be  analytically  simplified.  Fig.  2- 
2(b)  shows  an  equivalent  cavity  which  produces  identical 


When  the  effective  cavity  length  exactly  equals  one  cavity 
wavelength,  all  the  incident  energy  will  be  reflected  in  the 
-6  direction.  This  total  reflection  is  shown  in  Fig.  2-14. 
All  rays  entering  the  cavity  above  this  boundary  ray  will  be 
reflected . 


Fig.  2-14.  Total  Reflection 
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Summary 

There  are  several  steps  necessary  to  conduct  this  ray 
tracing  analysis. 

Step  1:  Calculate  the  peak  RCS  of  the  cavity  using 
the  following  equation: 

a  (peak)  =  4_tt_a2 

2 

X 

o 

Then,  for  each  incident  angle  0,  the  following  steps  must  be 
performed : 

Step  2:  Multiply  this  peak  RCS  by  cos^9  due  to  the 
varying  effective  aperture  area.  This  modified  RCS 
would  be  a  final  answer  if  all  the  energy  was 
backscattered. 

Step  3:  Calculate  the  cavity  wavelength  according 
to  the  equation: 


x  =  2w 
tan  0 

Step  4:  Subtract  the  maximum  integer  number  of 
cavity  wavelengths  from  the  cavity's  true  length  to 
yield  the  cavity's  effective  length. 

Le  =  L  -  ( n  x  )  rp.ax 


M>"  -H  'y 


btep  b:  Depending  on  the  effective 

calculate  the  percentage  of  backscatter 
to  the  following  equations. 


1  e 

acci _ 


L 


e 


%  backscatter 


0  <  Lp  <  X 

T 


<  LP  <  X 
7 


<  Lp  <  3  A 
T~ 


3_X<  L  <  X 
4 


2  Le  tan  9 
w 


2 


1 


Le  tan  9 


w 


2 


(Le  -  X/2)  tan  9 
w 


(Le  -  X/2)  tan  9) 
w 


Step  6: 


Finally,  calculate  the  RCS  as  follows: 

°  (©)  =  (  a  pea^)  (cos ^  9)  (%  backscatter) 
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III.  RCS  Measurement  of  a  Two-Dimensional  Cavit 


Introduction 

It  is  the  purpose  of  this  thesis  not  only  to  develop  the 
ray  tracing  theory  as  outlined  in  Chapter  II  but  to  formulate 
a  computer  program  to  accurately  predict  the  RCS  of  a  two 
dimensional  cavity  as  well.  However,  unless  the  computer 
predictions  can  accurately  compare  to  measured  results,  the 
theory  and  computer  program  are  without  value.  Therefore,  an 
important  part  of  this  project  was  the  measurement  of  the  RCS 
of  an  actual  two  dimensional  cavity.  This  data  was  then  used 
as  a  benchmark  against  which  the  computer  predictions  were 
compared . 

Design  of  the  Cavity 

The  primary  difficulty  associated  with  designing  a  two- 
dimensional  cavity  is  that  the  surrounding  structure's  RCS 
must  be  low  compared  to  the  cavity's  RCS  for  a  wide  range  of 
incident  angles.  One  such  design,  and  the  one  which  was 
eventually  used,  is  shown  in  Fig.  3-1. 

Embedding  the  cavity  in  a  wedge-cylinder  as  shown 
provides  several  advantages.  First,  none  of  the  supporting 
structure  is  normal  to  the  radar  except  at  large  incident 
angles.  And  second,  the  cylinder  behind  the  cavity  prevents 


a  large  diffraction  component  from  originating  at  the  back 


v> 


weight  could  not  exceed  75  pounds  due  to  limitations  of  the 
radar  range's  target  pedestal.  The  finished  aluminum  test 
cavity  weighed  an  unwieldy  66  pounds  and  was  constructed  from 
a  solid  cylinder  and  two  flat  plates  which  were  planed  down 
to  form  the  wedge  shaped  cavity  walls. 

All  the  seams  were  sealed  and  smoothed  to  allow  the 
propagation  of  traveling  waves,  and  the  two  leading  edges 
were  sharpened  to  minimize  the  two  edge  returns  for  vertical 
polarization. 

Test  Facility 

The  measurement  of  the  cavity's  RCS  was  performed  at  the 
Passive  Countermeasures  Branch  of  the  Air  Force  Wright 
Aeronautical  Laboratory's  Electronic  Warfare  Division  (Eldg. 
821)  under  the  supervision  of  Dr.  Robert  Puskar  and  Mr.  Ed 
Williams.  The  fully  instrumented  radar  cross  section  range 
consists  of  an  anechoic  chamber,  a  rotating  target  pedestal, 
a  computer  control  console,  and  several  different  platforms 
supporting  radar  transmitting  and  receiving  equipment.  The 
measurement  procedure  is  started  with  a  nulling  out  of  the 
return  of  the  range  itself.  These  range  returns  are  caused 
by  the  pedestal,  walls,  absorbing  cones,  and  other 
structures.  By  adjusting  the  magnitude  and  phase  of  the 
signal  (without  a  target)  in  a  second  receive  channel  which 
is  added  to  the  primary  channel's  signal,  the  range  return 
can  be  nulled  out.  Next,  the  power  level  of  the  received 


signal  must  be  calibrated.  By  placing  a  sphere  of  known  RCS 


IV.  Comparison  of  Results 


Introduction 

A  computer  program  was  written  so  that  RCS  values  for 
the  test  cavity  could  be  predicted  for  a  wide  range  of 
incident  angles.  The  computer  code,  along  with  complete 
documentation  and  details  of  its  operation,  is  contained  in 
Appendix  A.  The  results  from  this  program  were  compared  to 
the  actual  test  measurements  so  that  a  qualitative  judgement 
of  the  usefulness  of  the  ray  tracing  theory  could  be  made. 

Measurements 

Measurements  were  made  at  a  frequency  17.2  GHz  for  both 
horizontal  and  vertical  polarizations  because  this  was  the 
highest  frequency  for  which  the  test  range  had  equipment. 
These  results  are  shown  in  Fig.  4-1  and  4-2.  There  is 
obviously  a  significant  difference  between  the  two  graphs  for 
incident  angles  between  -15  and  +15  degrees.  This  could  be 
caused  by  several  factors.  First,  the  cavity's  two  leading 
edges  create  more  backscatter  for  vertical  polarization 
because  of  currents  which  are  set  up  along  these  edges.  And 
second,  horizontally  polarized  energy  causes  travelling  waves 
to  be  set  up  on  the  cavity's  sidewalls.  This  discrepancy 


highlights  the  fact  that  this  test  was  not  of  a  truly  high 
frequency  nature.  Ideally,  if  the  wavelength  of  the 
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electromagnetic  energy  is  small  enough  compared  to  the 
cavity's  dimensions  (particularly  the  width  of  the  opening), 
then  the  measurements  should  be  the  same  for  both 
polarizations . 

The  measurements  were  plotted  using  a  5  degree  averaging 
window.  Thus,  the  RCS  value  for  each  incident  angle  is 
actually  the  average  of  RCS  values  for  five  degrees  centered 
at.  that  angle.  This  averaging  technique  smoothes  out  the 
curves  sc  that  RCS  variations  don't  obscure  the  general 
trends . 

It  should  be  noted  that  for  angles  in  excess  of 
approximately  45  degrees  off  axis,  the  RCS  of  the  cavity  is 
dominated  by  the  return  from  the  large  flat  plate  w'hich 
serves  as  the  outside  wall  of  the  wedge  cylinder.  Since  only 
the  return  from  the  interior  of  the  cavity  is  of  interest  for 
this  thesis,  comparisons  of  data  beyond  these  angles  is  not 
made . 

Computer  Results 

The  results  from  the  computer  programi  are  plotted  in 
Fig.  4-3.  Since  the  entire  ray  tracing  approach  assumes  very 
high  frequencies,  these  RCS  values  are  independent  of 
polarization. 

RCS  values  were  calculated  for  incident  angles  at  one 
degree  intervals.  So  that  a  fair  comparison  to  the  measured 
data  could  be  made,  a  subroutine  was  devoted  to  generating  a 
5  degree  averaging  window. 
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Also,  calculations  of  predicted  RCS  values  were  only 
made  for  incident  angles  of  no  more  than  45  degrees  off-axis. 
The  ray  tracing  approach  neglects  the  RCS  contribution  of  the 
outside  of  the  wedge  cylinder. 

Comparison 

Fig.  4-4  shows  the  computer  data  along  with  the 
measurements  for  horizontal  polarization.  It  is  clear  that 
predictions  for  angles  of  less  than  ten  degrees  off  axis  are 
not  accurate.  Indeed,  the  shapes  of  the  curves  for  these 
angles  are  not  similar.  However,  beyond  these  angles,  the 
ray  tracing  theory  is  more  successful.  It  accurately 
predicts  peaks  at  angles  of  approximately  +  22  degrees  and  +  33 
degrees.  Nevertheless,  although  the  trends  are  the  same  for 
these  angles,  the  magnitude  of  the  RCS  peaks  for  the  measured 
data  falls  off  faster  than  the  cos2&  decrease  predicted  by 
the  ray  tracing  theory.  This  could  be  caused  by  the  fact 
that  as  the  incident  angle  increases,  the  effective  width  of 
the  cavity  decreases  causing  the  test  to  be  less  and  less  of 
a  high  frequency  case.  This  causes  the  theory  to  be 
decreasingly  valid  as  the  incident  ancle  increases. 

Fig.  4-5  shows  the  same  computer  data  compared  to  the 
measurements  for  vertical  polarization.  Here  the  agreement 
is  much  rr.ore  striking.  The  overall  shapes  of  the  two  curves 
are  very  similar  with  the  computer  accurately  predicting  the 
angles  at  which  RCS  peaks  occur.  The  agreement  is  also  quite 
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good  for  small  angles  in  contrast  to  the  horizontal 
polarization  case.  Once  again,  however,  the  computer  code 
predictions  do  not  fall  off  as  fast  as  the  measured  data. 
Nevertheless,  these  results  are  very  encouraging.  Evidently 
the  horizontally  polarized  wave  is  more  sensitive  to  the 
ratio  between  wavelength  and  cavity  width.  In  other  words, 
since  the  theory  agrees  more  closely  with  the  vertical 
polarization  case,  perhaps  the  measured  results  for  vertical 
polarization  more  closely  approximate  the  truly  high 
frequency  case. 
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V.  Conclusions  and  Recommendations 


Conclusions 

Overall,  the  comparison  between  the  computer 
calculations  and  the  measurements  was  encouraging  enough  to 
show  that  the  ray  tracing  approach  shows  some  promise  for 
limited  application.  From  the  outset  it  was  obvious  that 
this  approach  is  limited  to  high  frequencies,  but  this  does 
not  mean  the  theory  is  without  application.  Because  it  is 
straightforward  enough  to  be  implemented  by  a  short  and 
streamlined  computer  program,  this  technique  can  provide 
considerable  savings  of  computer  execution  time  and  storage 
requirements  when  compared  to  some  of  the  more  widely 
applicable  methods  described  in  Chapter  I. 

The  discrepancies  between  the  measured  and  calculated 
RCS  values  for  vertical  polarization  as  outlined  in  Chapter 
IV  are  probably  caused  by  one  of  two  phenomena.  First,  as 
mentioned  previously,  it  could  be  caused  by  the  decreasing 
aperture  width  forcing  the  cavity  to  approach  the  low 
frequency  case.  If  this  is  indeed  the  cause  of  the  less  than 
perfect  data  agreement,  then  it  demonstrates  the  lower  limits 
of  the  model's  applicability.  And,  the  computer  code  in  its 
present  form  should  become  increasingly  accurate  for  higher 
frequencies.  Another  possibility  exists,  however.  The 
progressive  gap  between  the  data  could  be  caused  by  a  failure 
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to  account  for  diffraction  of  the  radar  wave  both  as  it 
enters  and  exits  the  cavity.  If  this  is  the  problem,  then 
the  discrepancy  will  persist  for  all  frequencies  thus 
representing  a  flaw  in  the  theory.  Nevertheless,  a 
diffraction  coefficient  subroutine  as  described  by 
Kouyoumjian  and  Pathak  (3)  could  be  easily  included  in  the 
computer  program  without  serious  execution  time  penalties. 

Recommendations 

Since  the  ray  tracing  approach  to  cavity  scattering 
shows  promise,  it  may  merit  some  additional  work  in  the 
future.  If  further  work  is  deemed  appropriate,  the  following 
items  should  be  accomplished.  First,  future  tests  should 
either  use  higher  frequencies  or  larger  cavity  dimensions  to 
ensure  that  low  frequency  phenomena  do  not  distort  the  data. 
Since  the  theory  is  a  high  frequency  method,  any  tests  should 
use  frequencies  which  are  "high  enough."  If  this  does  not 
solve  the  problem  of  the  less  than  desirable  data  agreement 
for  large  incident  angles,  then  a  diffraction  coefficient 
subroutine  should  be  included  to  find  out  if  this  indeed  is 
the  cause  of  the  problem.  Furthermore,  different  models 
should  be  tested  to  make  sure  the  theory  works  for  cavities 
of  different  widths  and  depths.  This  will  serve  to  validate 
this  approach  for  a  general  cavity  rather  than  for  just  one 
specific  shape.  Finally,  and  most  importantly,  the  ray 
tracing  theory  should  be  extended  from  two  into  three 
dimensions.  Whereas  two-dimensional  cavities  are  almost  non- 


existant  in  the  real  world,  three-dimensional  cavities  can  be 
found  in  prodigious  numbers.  Therefore,  the  ultimate 
application  of  this  technique  depends  upon  its  extension  to 
handle  three-dimensional  cavities.  While  this  task  is  not 
trivial,  it  is  certainly  more  painstaking  than  complicated. 
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APPENDIX  A 


100  REM  THIS  PROGRAM  WAS  WRITTEN  BY  R. SCOTT  DERING 

110  REM  IN  AUGUST  1984.  IT  RUNS  ON  A  TOSHIBA  T-100 

120  REM  PERSONAL  COMPUTER.  IT  CALCULATES  THE  RCS 

130  REM  OF  THE  TWO-DIMENSIONAL  CAVITY  DESCRIBED  IN 

140  REM  CHAPTER  III  USING  THE  RAY  TRACING  METHOD  AS 

150  REM  OUTLINED  IN  CHAPTER  II. 

160  REM 
170  REM 
180  REM 
190  CLEAR 

200  DEFSNG  A-M.O-Z 

210  DIM  MAG  < 50  >  ,  ARRAY  <  1 00 )  ,  RCS  (  1 00 ) 

220  REM 

230  REM  FIRST,  COMPUTE  THE  PEAK  RCS  OF  THE  CAVITY 
240  REM  WHICH  IS  ACTUALLY  THE  RCS  OF  THE  SQUARE  FLAT 

250  REM  PLATE  AT  THE  END  OF  THE  CAVITY 

260  REM 

270  REM  LAMBDA  IS  FREE  SPACE  WAVELENGTH 
280  REM 

290  REM  AREA  IS  THE  AREA  OF  THE  CAVITY  END 
300  REM 

310  LAMBDA*. 02 
320  AREA* . 044 
330  W*.  1 
340  LENGTH*. 385 

350  F'EAK=4*3.  14159*  (AREA  2)  /  (LAMBDA ""2) 

360  REM 

370  REM  NOW  THE  PEAK  RCS  MUST  BE  SCALED  ACORDING  TO 
330  REM  THE  THEORY  OUTLINED  IN  CHAPTER  II. 

390  REM 

400  FOR  N= 1  TO  50  STEP  1 
410  REM 

420  REM  FIND  THE  "CAVITY  WAVELENGTH" 

430  REM 

440  THETA  -N/57.3 

450  WAVEL=2*W/ (TAN (THETA) ) 

460  L=LENGTH 
470  REM 

430  REM  DECIDE  WHICH  CASE  APPLIES  <1,2,3,  OR  4) 

490  REM 

500  IF  L  =. 25*WAVEL  GOTO  960 
510  IF  L  =.5*WAVEL  GOTO  990 
520  IF  L  =.75*WAVEL  GOTO  1020 
530  IF  L  =WAVEL  GOTO  1060 
540  L=L-WAVEL 
550  GOTO  500 
560  REM 


5IS  £13  5CALE  THE  RCS  SY  THE  CDSINE  SQUARED  FACTOR 
590  PERCENT  =F'ERCENT K  (COS  (THETA)  /  2 
aOO  REM 

610  REM  MAO ( N )  IS  THE  MAGNITUDE  OF  THE  RCS  FOR  ANGLE  N 
o20  REM 

630  MAG  <  N ) =FEAK* PERCENT 
o40  NEXT 
650  MAG ( 0 ) =PEAK 
660  REM 

670  REM  ASSIGN  RCS  VALUES  TO  NEGATIVE  ANGLES  AS  WELL 
630  REM 

690  FOR  N=0  TO  50  STEP  1 
700  ARRAY (N+50) =MAG (N) 

710  ARRAY (50-N) =MAG <N) 

720  NEXT 
730  REM 
740  REM 

750  LPRINT  "ANGLE  (DEG.)",  "RCS" 

760  REM 

770  FOR  N=2  TO  98  STEP  1 

730  LPRINT  " - " 

790  REM 
800  REM 

810  REM  AVG5  IS  THE  AVERAGE. RCS  OVER  THE  5  DEGREE  WINDOW 
820  REM 

•  330  AVG5=  <  ARRAY ( N-2 ) + ARRAY ( N- 1 ) + ARRAY ( N ) +ARRAY ( N+ 1 ) +ARRAY ( N+2 ) ) / 5 

840  REM 

350  REM  RCS(N)  IS  THE  AVERAGE  RCS  FOR  ANGLE  N  IN  DECIBELS 
360  REM 

870  RCS  <N) =10*L0G (AVG5) *.43429 
880  LPRINT  N— 50, RCS (N) 

390  NEXT 
900  END 
910  REM 

920  REM  THE  FOLLOWING  SECTIONS  COMPUTE  THE  SCALING 
930  REM  FACTORS  FOR  EACH  OF  THE  FOUR  CASES  ACCORDING 
940  REM  TO  THE  EQUATIONS  LISTED  ON  PAGE  36. 

950  REM 

960  REM  CASE  1 

970  PERCENT  =2*L* (TAN (THETA) ) /W 
980  GOTO  590 
990  REM  CASE  2 

1000  PERCENT  =2* ( 1-L*TAN (THETA) / W) 

1010  GOTO  590 

1020  REM  CASE  3 

1030  LF'RIME  =L-.5*WAVEL 

1040  PERCENT=2*LPRIME*TAN (THETA) /W 

1050  GOTO  590 

1060  REM  CASE  4 

1070  LPRIME=L~. 5*WAVEL 

1080  PERCENT=2*(1-LPRIME*TAN (THETA) /W> 

1090  GOTO  590 
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